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ABSTRACT 


The present work deals with the response of a 
cantilever pipe conveying fluid under transient co’npressive 
follower load The internal damping, external damping 
fluid flow, Coriolis component due to fluid flow and exter- 
nal compressive tangential load at the free end of the 
pipe are taken into account for the formulation of problem 
The partial differential equation describing the motion is 
reduced to a set of coupled linear ordinary differential 
equations with the aid of Galerlcln’s technique and the 
resulting equations are numerically solved by Haming's 
method 

The external damping stabilizes whereas the 
internal damping destabilizes the nonconservative system* 

The fluid damping and coriolis force due to it, overcome 
the destabilizing effect of internal damping* Whether or 
not the maximum response will increase or decrease as the 
load duration increases seems to depend on the deflection 
shape of idle column at the instant the load is released. 

For a particular load and load duration the second mode 
dominates the overall response of the column This indicates 
that in the case pf very high loads acting for a short 
duration the second mod© response also is of importance# 

(Phe magnification factor is a good measure of the overall 
response but does not always correspond to maximum 
response* 



CWTER 1 


INTRODUCTION 

The study of the response and stability of 
structural elements under nonconservative force has become 
increasingly important due bo the application in a variety 
of fit-lds After classifying the various types of forces, 
a brief survev of tie etistuig literature on the stability 
of nonconsorvative system is presented Only those litera- 
ture w'ich is relevant to the present vork is touched 
upon, 

1 ,1 Classification of Forces 

Forces vhlch act on any structure can be class!" 
fled broadly into active and reactive forces, Active 
forces c«an further be clissified is ohown In Fig 1 


Active^ I'oroes 

Non stationary loads 

5' (q, q, t) _ _ 

Velocity dependent 

^ (q> q) 


stationary loads 

1* <q» q) 


Gyroscopic 
dw = 0 




’ de'^endent Velocity independent 

q> q) Cq) 

non-circulatory circulatory 

1 dvr = " dv dw dv 

dissipative 

dw < 0 


FIGUBB 1, GLASSIFIC*tTIOH OF FQRG3SS 
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Nonst^tionaiy loads depend on time explicitly 
where as stationary loads do not The corlolis force 
due to fluid flow is thK, example of gyroscopic force, and 
material damping, external damping and fluid friction are 
the oxamplo of dissipative forces Under the category of 
velocity independent forces we have the noncircuLatory and 
circulatory forces Constant loads are the example of 
noncirculatory forces 

This classification enables us to put the stabi- 
lity problems of nonconservativn system, in throe distinct 
groups * The first is concerned with clastic systems sub- 
jected to the so called follower forces, that is, forces 
which follow in a prescribed manner tne deformation of the 
system Such forces called the “circulatory** forces, are 
nondissipative and in general nonconservative, A second 
group deals with stability of rotating shafts (whirling) 
and the third and last group with stability of elastic 
bodies placed in moving fluid (aeroelasticity) , In this 
work we are concerned with the first type of stability 
problem, 

1 .2 Destablllaing E»ffeot of Internal Damping 

Since the pioneering studies of Ziegler in 195rs 
on the destabilizing effect of internal damping of non- 
conservative elastic systems, many researchers have 
investigated the various aspects of the Influence of 
damping , 
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Neraat - Nasser and Herrmann [11, Wemat - 
Wass'^r [ 2] , Prasad and Hermann (3j and Bolotin and 
Zhl^ngex [ 4] have shown that the crlLical load of an 
undamped discrete or continuous, elastic system under 
follower forces is an upper bound for the critical load of 
the system with slight internal damping Jong [ 20] showed 
for the case of two degree of freedom systems that the 
flutter bound ^^y curves of the uiidampod systems are the 
envelope of the family of flutter boundaiT curves of systems 
with small internal damping The destabilizing effect was 
further elaborated by Bolotin , who considered a 
general system with two degree of freedom not related to 
any particular mechanical model and who found, additionally, 
that the des tabilizing effect 1 1 the presence of slight and 
vanishing dsunping is highly dependent on the relative 
magnitude of damping coefficients in the two degrees of 
freedom 

Additional insight iito destabilizing effects of 
linear velocity dependent damping in nonconservative system 
was supplied by < Herrmann and Jong[ 6 ] They stu- 
died the roots of Ui characteristic equation, in addition 
to the stability criteria and introauced the concept of 
degree of instability* The'^'’ established a relationship 
between critical loading for no damping, for slight damp- 
ing £«id for vanishing damping It was found that while 
the presence of small damping may have a destabilizing 
effect, prpper interpretation of tJio limiting process of 
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vanishing damping leads to the suae critical load as for 
no ds-iiping Herrminn and long [ ?] showed further that for 
certain systems instability b-^ divergence or by flutter 
may occur depending upon the ratio of the damping coeffi- 
cients in liie degrees of freedom 

however, the destabilizing effect of viscous 
damping in a linear, dynamic sysxiom with two degrees of 
freedom, subjected to nonconservative forces has nob been 
put in a strong theoretical found tion Nor it has been 
shown whether a more general svst-ra with many degrees of 
freedom can also exhibit such beiaviour* Moreover from 
examples v/orked out in [ 5 , 6, 7l xt is not clear whether 
other velocity dependent forces t8,9l can have similar 
effects I 

Nemat - Naseer and HerrraEinn[ 10 1 have proved 
that in a genc.nl nonconservativo system vdth N - degrees 
of freedom not only slight damping but all sufficiently 
small velocity'' dependent forces imy have a destabilizing 
effect 

This oifect in continuous systems was considered 
hy ' Nemat - Na-sser, > Prasad and Nerrmann C 111 
The velocits^ dependent forces raiy be due to either internal 
or externa] dajnping, or to corlolis forces (in pipes 
cohveying fluid) or to other gyroscopic forces* For 
continuous systems with slight damping, Nemat - ilasser 
ft §1 proved that the flutter load parameter of the un- 
damped system is an upper bound for that system 
vlth slight damping, 



1 3 Effect of Extprnal Dajnping 




file effect of external damping on a linear elastic 

cantilevered column subjrcted to a constant follower force 

at lbs fro^ iid, is not des babilizing , has heen shown by 

Nemat - lass^r, Prasad and Herrmann [ 11 ] It is s’lown by 

Plaut and Infante [12 1 that the critical load increases \d.th 

"W T 

increasing dampiUi froin vlie value of 20 05 -n at zero 

1 *^ 

damping to th» limiting value of 37 7 ^- 5 - for large damp- 

ir 

ing This behaviour is in marked contrast to th-it of the 
internally damped column and aljo to that of the conserva- 
tive systems (where external daiping has no effect on the 
critical load) 


1 4 Stability of Pipes Containing Plowing Fluid 

lienzamin ( I 3 1 has investigated that, when ttio 
tube is fixed at one end and free at the other, large flow 
velocities can cause instability in certain modes of vibra- 
tion, whoreas for smaller velocities all possible modes are 
damped* 

tow IS s\«nptions were given by Benzamin for the 
formul‘=ition of problem of solid and fluid interaction If 
the lengths oi the pipes are much larger 'Uian the dimen- 
sions of t’ exr cross-sections, it Is reasonable to assiame 
that the small scale details of flow, for example, turbu- 
lence end secondary flows arising at the bonds in the 
flow passage will have an insignificant effect compared 
with that of the mean flow* Hence tne forces exerted by 
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the fluid are of two kinds First, fcliere are transverse 
forces which arise from displacements, this is in effect 
daa to a linear momentum flux bangonbial to the pipes 
These forces 1 nclude the coriolis reactions and reactions 
due to transverse acceleration of the fluid mass Secondly, 
there are forces due to fluid friction B\is category 
must include both the taigoiitial she^r forces and normal 
forces due to bi(» pxesjure developed in the fluid as a 
result ol friction The important conclusion is that these 
forces have a neutral octal effect on the pipr** Thus the 
dynamical problem is independent of fluid friction 

An exhaustive literature survey vras given by 
Gregoiy and Paldoassis ^ 141 and PaSdoussia 
and Issidt 15 * 1 from w’lich we see a remarkable develop- 
ment of the subject in tlie twenty years sjnee the Trans - 
Arabian pipeline was observed to vibrate presumably as 
a result of Interaal flow However, the impetus for most 
of those studies did nob come from a desire to solve a 
practical problem, but rather from one or more of the 
following considerations 

(i) the physical problem was inher'’ atlv intriguing, 

(ii) the problem offered scope for •’nteresting mathematical 
manipulation, (iil) the problem represents one of tiie few 
oases where a Cnon-oonservativo) follower load is jhysically 
realizable f The force of this last poin-t become^ obvious 
upon recalling bhe very considerable tlieoretical interest 
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in the stability of column subjected to a tangential 
follower lon.d ao the free end [ 16 ] ^ even when no means 
was ovidoat for producing it N'ow, internal flow existing 
from the free end of a cantilever docs generate a tangen- 
tial load Tnis being a physicall'''’ easily realizable 
system - as compared to usung solid - fuel rockets Quite 
recently, however, increasing iiidustrial iiiterest has 
developed in the dynamics ol pipes convt.-»?ing fluid, in 
conjunction with detormining the rosponae of such systems 
to arbitrary force fielas, with api^lication to heat 
exchanger, liquid- fuel rocket piping and nuclear reactor 
coolant channels 

1 5 Effect of TrinsienL Follower Load 

S«T Nouh and V bundararaoan I 17lhavt- studied 
the ampllf ic ition of the lateral response, ol a cantilever 
column assuming a sot of ijiitinl pf^rturbations, effect of 
pulse duration, internal damping and the effect of diffe- 
rent init1 ’1 perturbations I hafvo sasso been dise^ns-sed. 

1 ,6 Present Work 

fhe response of ^ linear olastlc cantilevered 
system subjected to a lonconservative load has been dealt 
with in this thesis The following spr'cific problems 
hav6 besn studied 

1 . Cantilever column subjected to a tangential follower 
load The affect of exte^hal and internal damplhg 
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on the first critical follower load have been studied 
The combinf d effect of both dampings on critical load 
is ilso discussed 

2 OaJitilover pipo conveying fluid under tangential 
follower load The elf act of Iluid flow, coriolis 
lorce du j to lu and internal d-iiping on oxtornal 
tangenLial Aollower load (loS'a than idle first critical 
load 20 05) have heei discussed 

3 Oantilover column under triaisiont lollower load 
The effect of very high loads (several times the 
first critical load), acting for a short duration, on 
the response of the column is discussed The effect 
of load duration on the response of column is also 
studied* 

Cantilever pipe conveying fluid under transient 
tangential follower load The effect of fluid flow, 
coriolis force due to it and internal damping on the 
response of the pipe subjected to high loads (several 
times the first critical load), is studied 

Tho equation of motion is found by equilibrium 
method In the present work the stability problem is 
discussed from tho response of thu cantilever column* The 
Gilerkin's technique is used for reducing the linear 
partial differential equation to ordinary linear diffe- 
rential equation and Homing's method is used for perform- 
ing the solutiort of coupled linear differential equations 



CHAPTER 2 


FORMULATION 

In tills chapter the equation oi motion of a 
cantilever pipe conveying flu d under external follower 
load is formulated by the equilibrium method 

2 1 LquaLioa of i otLon 

A uniform cantilever pipe of length L, flexural 
rigidity El and mass per unit length mp, conveying a stream 
of incorapres ible fluid of mass m^, per unit length with a 
mean flow velocity U and a follower force Po acting at 
its free end is shown in Pig 2 (a) The x -* axis coincides 
with the centre line of the undeflected pipe md its trans- 
verse deflection y is measured perpendicular to ox Both 
ox and oy lie in a horizontal plane 

The free body diagrims of fluid element and a pipe 
element are aiown in Fig 2 (b - c) 

The cross sectional area through which the fluid 
passes is A, the iitemal perimeter is S and the internal 
pressure is p Consider now a small element of the pipe, of 
length fix, and the corresponding element of tiie enclosed 
fluid, of volume 5D The rate of change of momentum over 
fi B may be written as 

^ s Jj7 [ If + ( ) V ] P dD (2 1) 

Where ^ is the instantaneous velocity of a small element 
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dD within 6 D Assuming that the radial variations in the 
flow velocity are small and that secondary flow effects are 
negligible - which was actually presupposed in "ttie assump- 
tion of planar motions iaiat is, assuming that one lias 
plug flov/ and hence ^ may be written as 

5 +3-^ = tj i +1^ 5 (2 2) 

Moreover, upon assuming that y and are small 
and that motions are of long wavelength such that is 
small, Fquations (2 1) and (2.2) after some manipulation 
yield 

f|=®fH ^ U |-x ^ 

(2.3) 

Hence, for the fluid element of Fig 2 (b), force balances 
in the x and y directions yield 

A + <lg S = 0 (Z.k) 

y = 0 (a.J) 


where q- is the shear stress on the internal surface of the 
pipe and F the transverse force per unit length between 
pipe wall and fluid. 


The equations of motion for the pipe may similarly 
be found as 


is 

dx 




Q 


3^ y 

3x2 


0 


(2.6) 
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and 




dX 


+ T 




0 

(2 7 ) 


where is the external follower load acting at the free 

end of the cantilever, T is the longitudinal tension, Q is 
the transverse shear force in the tube and K is the 
coefficient of viscous damping due to friction of iAie pipe 
with the surrounding stationary fluid medium The shear 
force Q is related to bending moment acting on tJie 
section of the tube by 


a 

+ Q = 0 (2 8) 

provided rotary inertia is neglected, while is related 
to the lateral deflection by 


M. 


= E I 


3 X 


From Equations (2*8) and (2 9)> we have 


(2 9 ) 


5 -ni ( 210 ) 

Alt 

Wo consider B as the coefficient of internal 
dissipation waich was assumed to be viscoelastic and of the 
Kelvin - Voigt tyoe Then 


Q _ - ( -^'+ E) I (2.11) 

The terms of second order magnitude are neglected- 
in accordance with the Buler - Bernou31i beam approximation 
for small lateral motions On combining B qua tions (2 
(2.7) and (2.11 ), we have 
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9^ y 9 


<2 12 ) 


q. may “be eliminated between Eq (2 and Eq, (2 6) 
s 


to give 


(p A - T) =: 0 


(2.13) 


Thus (p Ji. - T) is independent of x, furthermore, 
at the free end p A - T = 0 and consequently 


p A - T = 0 


(2 1 ^) 


everywhere * This is an important conclusion since it 
establishes •Uiat the dynamical problem is independent of 
fluid friction. 

On substituting Eq (2 1^) into Eq» (2.12), we 


obtain 


( + E) I 2^ + mj, 4. 0 7 




(2 1 ?) 


"'’i g* ^ g I (P *02) ^ 

&x^9 t 9x^ ° 


0 


( 2 . 16 ) 
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This is tho required equation for studying the 
response of a cantilever pipe conveying fluid under follower 
force with tho boundary conditions as follows 


and 


y = 


a x' 


-.a-Y, s: 0 
a X 


at X = 0 


3 

= ^ 
8 x*^ 


(2 17 ) 


= 0 at X = L 


2.2 Nondi niensionalization of the Equation of Motion 

Before we proceed with the analysis, it is 
desirable to express the problem in dim' nsionless parameters 
and accordingly we choose 


^ = 5 . ^ 5 » 


= ( 


B I 


“f 


) 


1/2 


2 » 


. Mf. 1/2 

u =: ( g±^) UL 


r .2 - 

, r = 


m^. 

+" 


_ E / I K 

" r 2 ^ E W + ra^9 * ° » 

^ ^ P 2(B I (mj. +mp)) 


1/2 


and 


P. L 


p ( T ) =: >0 in CO , T) and P ( t) = 0 when 

ii ^ 

T >■ T , whore ^ is the load duration 

Now Equation (2 16) can be written in terms of 
dimensionless paramo ter d as 


3 ^ w . _3^ w. 

" L. + * 


3? Vt 


+ {p(t)+u^}^ ^ 


3 w 


+ 2 u r +Zg . 

elCSf 


^ 3^ w 
‘'3t2 


3€ 


= 0 


(2,18) 
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and. 'boiindaiy conditions reduce to 


32 ^ 




T? 


at 5 = 0 


at 5 = 1 


(2 19) 



CHAPTER 3 


SOLUTION 


Tho solution procedure for Eq (2.18) is discussed 
in this chapter A brief computational approach for solving 
a set of ordinary coupled differential equations is also 
given 


3.1 Solution Procedure 


Equations (2 18) arc reducod to an infinite system 
oi ordinary differential equations through a Galerkin's 
procedure We now apply the Galerkin's method for the 


analysis of system (2 18) 
eigenfunctions { ( ^ ) ) 

ing eigenvalue problem 

■? 


Wo consider a set of orthonoimal 
, obtained by solving the follow- 


A 

d <|> 


m 


.il ^m 


= 0 


^ ‘I'tn 

*m ~ ar” = 


at 


d^ it 


m 


d^ «)> 


m 

T 


5=0 

0 at 5 si 


(3 1) 


dS"^ d 5' 

The orthonormal eigenfunction or characteristic 
function obtained from ’’’‘quation (3 1 ) is given below 


( 5 ) a { Cosh ( 5 ) - Cos (3^? ) ^ 

[ Sirih ( Sjjj 5 ) - Sin ( )) } (3 S) 

Cos ( 0 ) + Cosh (6-,) p 2 
where a„ = % 
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and satisfies the frequency or characteristic equation 

Cosh ( Cos ( Bjjj) =5 -1 

The values oi p and a are gxven in a tabular 
form in ^puendix (A) for various values of m 

Let us assume a solution of the form 


w (C,t) = I 4> (5) 

n=1 ^ ” 

and variation in w of the form 


CO 

® ^ ~ m=1 

for Eq (2 18). 


(3 3) 


(3 


Substituting Eq. (3 3) into Eq (2 18) and multi- 
plying botSh sides by 6w, we get 

« d^<|> o " d^<jL 

1 I?q„ (t)6w + {P (t)+u^} I (t)6w 

n=1 d r ^ n=:1 dr ^ 

00 , 00 
n=1 ” ^ n=1 d ^ 

+ 2 u r I a (t) 6w+2c I (C)q„ (t)6 

= 0 (3*5) 

where a dot Indicates the differentiation witti respect to f 


w 


Substituting Eq. (3 into Eq (3 5) and integrating from 
25ero to unity we get 
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1 » 

d C 

00 

2 


/ X 9" , n ( J X d) (5 ) 6 a ( T ) 
1 n =1 j t-M- ' wsl m ^ ' 


d5 


d% 


+ / { P ( T ) + u^) X “3 Si ^ ^ ^ ^ ^ 

« n=1 d 5 m=1 


d? 


o 

+ s' i ♦„ (5 ) (L, (t ) ! ♦_ «(J- (T ) ac 

n =1 ^ " ni =1 “ ^ 

0 

1 , «» dN> *** 

+ / V X jpr % \ V ^ ®Sn ^ 

Q n=1 d 5 in=1 

1 “ d<I) * 

* s a " -xf '»n ' ^ ^ Jl ^m ^ 5 ’ « 

o 


+ 2 0 / J ♦n ^ a_ (f ) (5 ) « %, 't ) -If = 0 

Q n =1 m =:1 


or 


\ + L ‘mn + P < ^ ♦mn + ♦jin^ % 


00 00 . ^ 
nil ^ nil ^ 

- a 0 i U <>„ = ° 

where i> 


<l> 

^ mn 

1 

- f 


(C ) 


•d5 

^ mn ^ ^nm 

^ f 

0 

‘^m 

d 5^ 

mn 

1 

= / 

^ nt 

^ in 

(C) 

d<fn 

dC 

dS 

Ip il: ^ 

^mn ^ nm 


0 






6 

mn 

1 

- ; 

^ in 

(C ) 

(5 ) dE a 

1 m 

for 

0 m 


(3 6) 

for m n 

for m / n 

= n 
^ n 
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Ultimately '"iq.'uation (3 6) reduces to 
q^ + {P(t)+u } 

’** ^ ^ n Kn% ^ m = 1, 2, 3, 

(3 7) 

where ^ are the roots of the frequency equatioi of the free 
undamped vibration of tne beam 

3*2 Computational Approach 

Haming's method (Predictor - corrector method) I 18 1 
has been used for performing the solution of coupled Jlqs* 

(3 *7)* Since this method solves a set of first order diffe- 
rential equations only, Bqs (37) are reduced to a set of 
first order differential equations The redaction procedure 
is given in Appendix (B) Equation (3 7) contains two 
integrals and which are computed numerically as 

well as emperically The e iporical formulae for finding 
the integrals and , are given in Appendix (0) 

The convergence of GalerAin's method for infinite 
system of Bqs» (3 7) is checked by increasing m and 
observing the asymptotic nature of tho solution The 
convergence of Haming's method is simultaneously checked 
by varying the integration step size and observing the 
asymptotic nature of the solution The accuracy of the 
solution is seen up to fifth place of decimals ^ The 
double precision has been used for performing the solution 
of coupled equations. 



CHAPTER 4 


RESULTS AND DISCUSSION 

The numerical solution of the coupled equations 
(3 7) is performed for tie following cases 

1 Effect of damping on critical follower load for 

a cantilever column 

2 Cantilever pipe conveying fluid with lollower load 

3 Cantilever column under transient follower load 

•f Cantilever pipe conveying gluid under transient 

follower load 

The results obtained in tlie form of dimensionless 
parameters are discussed in this ciapter The effect of 
various dimensionless parameters on the response and 
magnification factor in the above cases have been discussed 

Initial conditions of iaie following form are 

assumed 

(0) = 10“^, q^ (0) = 0, m > 1 (4 1 ) 

A«} a measure of the total response, the magnifi- 
cation facbor M is introduced, defined by 

1/2 1/2 

M = { r (t ) > / q^ (0) ^ (4 2) 

The convergence of tie GalerKin’a procedure for 
a continuous system is checl sd by considering different 
values of m. The value of m = 5 iS finally used in the 
solution of Bqa (3*7) coupled jSqs (3 7) are 

integrated numerically by Haming's method (Predictor - 
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Corrector Method) The convergence of Hairing *s method is 
checked hy observing the numerlcil values up to the fifth 
place of decimals for two successive integration step sizes 
The integration step size of 0 002 is finally used for 
performing the incegrataon of Bqs (3 7) in all the above 
cases 

4 1 nrfect ol Damping on Critical 
iollower Load (u =: 0, r = 0) 

It is well known that in the case of a cantilever 
column subjected to a tangential follower load, the critical 
value of the exLernal load is 20 0^ The corresponding 
critical load in the case of a conservative load is 2? 

The effect of external damping on the oritical follower 
load i s studied from the response behaviour of the column 
The critical load increases with increase in the external 
load, which is 38CFig 3) This result checks with tiiat 
of Plaut and Infante ri2]|w^o have established liiis 
considering the problem as an eigenvalue problem Thus 
the external damping stabilizes the system 

The effect of various values of internal damping 
on the critical load is shown in the Fig if. As internal 
damping increases the critical load decreases, and for a 
particular value of i the critical load reduces by 
almost 60/^* Thus the effect of internal damping is 
destabilizing, which is in contrast to that of external 



21 


damping For a given follower load, therefore, an increase 
in internal damping mav causo a stable column to become 
uns table 

Assuming the same value of internal and external 
dampings as 0 001 an attempt is made to explain the con- 
trary beiaviour of these two dampings Xhe magnification 
factor vrilh internal damping, extornai damping and without 
damping are given in Table ^ 1 Fy closely studying the 
magnification J actor it can bo noticed that the effect of 
internal damping becomes more domljiant after an elapse 
of time T as compared to that of external damping It 
seems as if ttie coupling due to associated with the 
internal damping is responsible for the unusual behaviour 
of internal damping on critical foJ lower load, whereas the 
external damping is not associated \rt.th any coupling of 
tne modes 

The combined of loci of both external and internal 
dampings on critical follower load is shown ia Table 4- 2 
With an internal damping of 0 01 the critical load decrea- 
ses from 20 0? to 11*00 in absence of external damping* 

V/ith a consi,aat value of Y = 0 01, the increase in external 
damping (0=1) increases the critical load from 11 00 
to 19*00* reason for this behaviour of damping lies 

in the relationship between the restoring forces caused 
by internal and external damping, both of which depend on 
the rate motioh but in different manner^^ Ail 

4Mb Nk> rv 62131 ^^ 



TABLUI h 1 
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Tho variaUon of jna£ nif icntion factor vjlth time (at 
sJightly loss i/ha,D th first oritio'^l load, 20) for 
internal damping, txtf jnal daraping 3 nd i/ithooit damping 
is givon bx low 


n = r = 0 


Timo , T 


0,^ 
0*8 
1 2 
1t6 

2»p 

2,4 

2t8 

3*2 
3 6 

4.0 
4 4 
4.8 
?,2 
5*6 
6 0 

U 

I" 

6.0 


[7=? 0.001 


2.046 

1.696 

2,772 

f’S 
11 ,087 
10,304 

1 633 
17.721 
23 ?33 

7.37? 
25,060 
47 114 
30 371 
27 301 
84 913 
81 797 
11 .077 
13? .365 


1°. 


H with 
o s 0*001 


H with 


2.079 

1 419 
3,892 
6 029 

0.996 

8,306 

4,386 

^*760 

7,5‘+4 

0 559 

?,9p2 

2.749 

Vtl 

1 136 
2.788 
0.990 
5,032 
5 714 

2 602 


..L ° ° 


2.079 

1 ,420 

3.897 

6 038 

0 99? 
8,326 
4 398 
5.780 

7 571 
0 562 
5,928 
2,762 

1 505 

1 355 
1 143 
2,806 

m 

5.757 

a.62<f 
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TABLE h 2 

Tht, oiifjcb of oxti rnal dnjuplng on critic il followtr load 
Vov^ping tco internal damping constant, is shov/n "bfelow 

u = r = 0 


IntJrnaT'' 

5 Extlrnal 

( Critical 

daHiplng 

i damping 

0 follower 
^5 load 

0 01 

0 0 

11 0 

0 01 

0 01 

11 5 

0«01 

0 10 

14 0 

0.01 

1 00 

19 0 



external damping may decrease tihe effect of internal 
damping and the total restoring force on the column may 
bG reduced 

If 2 Cantilever Pipe Conveying Fluid 
with Follower Load 

In this case the efiect of fluid flovir and inter- 
nal damping on external tEUigentlT,! follower load (less 
than the tirst critical load 20 05) have been discussed 
Tho load is assumed to be 15 00 The variation of the 
magnification factor with timt is shown ii Fig 5 for the 
cases with and without fluid flow and internal damping 
In the absence of fluid ilow (u 0) the internal damping 
destabilizes the system, which should be the case since 
the load is f^r ater than first critical load with Y = 0 01 
For different value of u the effect of flowing fluid on 
M is shown in i'ig, 5 Witn a fluid flovr of u = 1 it can 
be seeii that tlie system is damped dovni anu the magnifica- 
tion factor decreases rapldlv with time Thus the fluid 
damping and coroilis force due to it, overccmie the destabi- 
lizing effect of Internal damping and stabilize the 
system. 

Thcro exists a critical fluid velocity, for a 
particular value of the mass ratio and external load, for 
which the system becomes unstable For example when the 
mass ratio is 0.295 "the external load is 15 00 the 
orj tical flaid velocity is ^ 5 
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k- 3 Cantilovar Column Under Transient Follower 
Load (u = 0, Y = 0, c = 0) 

The effect of load duration on th.o response of 
a cantilever coDumn under a traisi-nt follower load of 
magnitude P = 27 is shown in Fig 6, 7 Since the 
response is dominated hy first mode, only tho response 
‘l-j ( ‘’^ ) / (O) is plotted against tine t The duration 

for which thu load acts is denoted by x It can he seen 
that the maximum amplification of th» initial disturbance 
depends on tho Limo at which the load is released and 
always occurs after the rplc ase of the load For a load 
duration from 0 9 to 0 975 the response decreases with 
increase in load duration* As we incroese the load from 
0 975 tho response starts increasing (Fig 7) 

The vnriition of the maximum magnification factor 
(M ) with the load duration x is suown in Fig 8 Whether 
or not the maximum response increases or decreases as the 
load duration increases scorns to depend on biie deflection 
shape of tho curve at the instant th? load is released* 

This cm be s^en from Figs 9 (a - d). 

The effect of the load duration on the magnifi- 
cation factor in the case of a constant external Impulse 
of magnitude 21 is shown in Fig 10 With the load 
duration x = 0 3j M reaches a value of tho order of 
10^« With r* s 0,5} M*** decreases to ^70 end with 

)|c 

T «s 1 to 11, Thus, the laximura magnification factor 
decreases considerably as the load duration increases 
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An interosting case has been plotted in Fig 11 
which corresponds to P = 120 and t' = 0 2, In tiis case 
it cm be soon that liie second mod3 dominates the overall 
response eventnough tiie Dnicial condition is assumed to 
contain the first morie onlv This ■cype of phenomenon may 
be expected in the case vfhen tie external load is very 
high (several times the first critic'll load) and the load 
duration is smll as compared to bh natural period of the 
column. The study of this behaviour is made by keeping 
the load duration constant (say = 0 2) and varying the 
external load It can be seen from the I‘igs 11 , 18 » 13 
and 14* that for values of P greater or eq.ual to 120 
(P > 120) ttic second mode dominates the overall response 
of the oolumn eventhough the initial condition is given in 
first mode only. However, for a particular load duration 
T there exists an uppi r limit of the load beyond which 
this phenomenon disappears This can be seen from the 
Fif. iH- which cor osnonds to the load P = 135. For a given 
load d\iration it scorns that there exists a range of load 
in which Lht second mode dominates the overall response of 
the column. At very higli loads one should be aware of 
the responses in higher modes also eventhough the initial 
perturbation is assumed to contain the first mode only* 
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4 4 Cantilever Pipe Conveying Fluid 
Under Transient Follower Load 

We discuss in tiiis case the response of a pipe 
conveying fluid under a transient tangential follower 
load greater th^i i first critical load acting for a short 
duration t Th*’ effect of internal damping, flowing 
fluid, coriolid force due to flowin^ fluid and combined 
effect of nil ticse have boon shown ii Fig The 

external load is assumed to be P = 70 and Lho load acta 
for a duration of t** = 0 3, The maxi»mara magnification 
factor (M ) when no internal dampin^ is present (y = 0) 
is seen to be verv high, which should be the case since 
the load is more tlian three times the first critical load 
M* in this case is of the order of 1 0^ With y = 01 , 

m'* decreases to a value of ^30 indicating the effect of 
internal damping in reducing M in the case of a transient 
follower load. Whereas in the case of a steady load 
(slightly less than first critical load without damping) 
the internal damping destabilizes the system ( that is, 
the magnification factor increases with time) and in case 
of a transient follower load many times the first critical 
load acting for a short duration the effect of internal 

>v 

damping is to raducs M ®ie effect of flowing fluid is 
to reduce M further This can be attributed to effect 
of the centrifugal and coriolis forces due to the fluid 
flow When both the internal damping as well as the 
flowing fluid are present, then M reduces considerably, 
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The effect of mass ratio r ( ^ — j- _) 

mass of fluid + mass of pipe 

on M is also shown in the Pig 1? The increase in mass 
ratio will causon.decraase in M 




It can he noticed that 
in tjie case of fluid damping^ M' occurs much quicker than 
bhat in the case of internal damping 


Figures 16 and 17 correspond to the case where 
the second mode dominates the overall response ot the pipe 
even tliough the initial condition is assuraed in bho iirst 
mode only The respons'** ( t ) / q^ (0) is plotted 
against t for differ nt modes with and wilSiout flowing 
fluid. The load is assuraed to ho P - 120 and load duration 
T s 0,17^ The presence of flowing fluid {u = 1 ) does 
not cnange thw maximum response in tic second mode The 
second mode still dominates tie overall response All the 
other modes (including tac first mode) die dovm rapidly 
compared witn the decay o± the second mode* 


1+,^ Comment on the Hagnification Factor 

The magnification factor was defined to represent 
tho overall response of the column and the moxim^om magnifi- 
oation factor gives an idea of tho maximum response But 
ill siould he pointed out tiat -Uio inaximum raagnificabion 
factor does not alwaj^s correspond bo th' maximum response 
Tl^ls can ho soon from Fig 9(c)« however in most of tho 
ca^es tho magnification factor is a good measure of the 
ovfirall response 



CHAPTER 5 


CONCLJSIONS 

Prom the study of Uio response of a lint-ar elastic 
cantiloverod system subjected bo an external follower load, 
the following conclusions are arrived at 

1 The internal dampinb destabilizes blio no aeons crv^tive 
systom The critical load docreasi^s uith increase in 
internal damping The external damping stabilizes the 
nonconsorvative system This b haviour is in marked con- 
trast to thao of the Internally damped column and also 
to that of the conservative system (wiore external damp- 
ing has no effect on the critical load) It s ems as if 
bhe coupling due to associated ^<lth the internal 
damping is responsible for the. unusual behaviour of 
internal damping on critical follower load, where as the 
external darapinc is not associated with any coupling of 
tho modes 

2, In case of external foilower load (slightlv less than 
first criLbici] load) the fluid damping and coidolis 
force duo bo it ovorcorao the destabilizing effect of 
in born al damping and stabiliz*-* the system There exists 
a critical fluid velocity for a particular value of the 
mass ratio and external load, for which the system 
becomes unstable. 
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3 It is found tnat tho aiaxiraura ampllfic ’tion of tt\e 
initial disturbance depends on thu time at wnich the 
load is r=>lcased md always occurs afuor tlie release of 
the load V/heJbher or not tht maxiarum response increases 
or decreases as tiv load duration increases seemed to 
depend on thr> deflection siapu of Uie column at the 
instant thf load is rol< ased In case of high loads 
(several tirnos llae first critical load) tho maximum 
magnilic ition factor decreases consider 3,bly vith increase 
in load duration 

4* An intoresting case has been established, in which for 
a particular load duration a ranj of load exists for 
which tiio second mode dominatr s tie overall response of 
bho column, tventhough tho initial perturbation is 
assumod bo contain the first node only Thus it can be 
concluded that in cast of vorv high loads (several times 
tho first critical load) acting for a short duration 
one should be awaru oi tho x espouses in hi^ei inodes also 
In this casu lie effect of fluid flow and coriolis forco 
duo to fluid flow is also studied Tho second mode stdll 
doniiiatos tho overall response Tho response in other 
modes (including first inode) damped out faster as 
compared to the decay in second mode response 

The fluid damping, coriolis force duo to it and material 
damping roduce- tho iiiaximum magniflc tion factor consi- 
derably as compared to the case with fluid or material 
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dampings alone In case of fluid damping M occurs 
quicker as compared to the case with the Internal 
dajnping 

6 The maximum magnific >tion factor does not always corres- 

I 

pond to the maximum response However, in most of the 
cases it is 1 good men sure of the overall response 

Future Vfork 

The unusual behaviour oi the internal damping on 
critical follower load should be studied rigorously 
Further thought should be given to this problem by consider- 
ing the loss factor, strain energy etc of the system In 
case of very liigh loads (several times the first critical 
load) acting for a short duration, the second mode dominates 
the overall response for a particular range of loads A 
theoretical explanation is needed for this special behaviour 
The transient external follower load of other types such ap 
triangular and blast type'idll be of Interest 
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The Values of and Involved in Equation (3 2) 


m 

1 

1 "in 

1 

1 .8^?^104-1 

0 734-0955 

2 

If. 69^0911 

1 0184-664.4- 

3 

7 85^7lf7^ 

0,999224-50 


10 9955‘+07 

1 ,000033553 

5 

1 Jf. 1371 684- 

0 9999985501 

6 

17.27875957 


7 

20 4-2035223 


8 

23*56194^88 


9 

26,70353753 


10 

29 84-513018 



Eor m > 5 

Since, p^ = (2 o -i 1 ) | “ 1 .0 
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Til© second order coupled linear ordinary differen- 
tial Eq.s* (3«7) are reduced to first order coupled linear 
ordinary differential equations for performing the solution 
by Hajning‘s method 

The second order differential equations are 
^ V «n 2 u r 

+ aoq„+sJjqjj = 0, m = 1,2, 3, 5 (B - 1) 

For m « 1, Bq* (B-1 ) can be written in expanded form as 


<1 + Y pI}’ q^ +{ P ( t ) + u^ > { q^ + '^ ^2 ‘la ■*■ *** 13 *^3 
+ 'J’Hi. ^ > + a u r q^ '**’*'12 *12 ^3 

+ ♦iv % +1’1J Ij }+2c + filf = 0 

Similarly other equations can be written in expanded form* 

Let q^ « y^ , <L^ = « y^, qii^ = y^, q^ - 7 ^ 

# • e f 

** ^2 ’ ^^2 ~ ^4^ *^3 ~ ^6* ^ * ^8’ *^5 “ ^10 

then Bq* (B f 1 ) can be written in the following form for 
different values of m, 


’ll 


d y^ 

dT 



y yi> 


• • •» 



^2 


CB - 2) 



54 


d 72 

= dr"' = ^2 yjn^ 



= -V e!^y2 

-{ P (x 

) + u^> { <^ ^ ^ y^ 

+ 7i2 73 



+ *13 75 + 


+ >>15 79 } 

- 2 u 

r{7i1 

yp 



+ *12 74 + 

■ ^3 ^6 

+ *14 78 -I 

^ ’*'15 yio> “^0 ya " 

• t>^ 71 







(B 

- 3 } 

«2 

- ^3 - f 

"* dT - ^3 

(y > yi> 

> 7„) 

= y4 


(B 

- M-) 

I 2 

~ dr “" ^4 

( f » yi » 

7ni^ 






=702 74 ■ 

( X ) + 

u® } f 7 2^ 

y^ +<^ 

22 y3 

+ <!)23 



+ *24 77 +*25 79 * 

- 2 u r{72i 72 

+♦2274 



+ *23^6 +'•’24 78 + 

* 25 7l c? 

-;c yi^ 

6 ^ 
“ ^2 

^3 








(B 

- 5) 

“13 

- - f 

(t , y^, 

> 7m^ 

= y6 


(B 

- 6) 

• • 

«3 

- ii6 ^ f 
” dPr ~ ^6 

( T , , 

• > yffl’ 






= -7 0^76' 

-{P (•»* 

) + u® }{ 7 

31 3^1 

■*■**’ 32 

ys 



+ <^3 75 + 

^^34 ^7 

+ 73579 } 

- 2 u 


I y2 

’"32 y4 


+ ♦aj 76 + ♦ 31 ^ yg +1’35 yio> ye - »3 ^5 

• ^ ^7 

\ = "dr = ^7 <^ > yi> 


> = ^8 


(B - 7) 
(B - 8) 



% 


= = ^*8 , y^, , y^) 

= “Y 8 4 yg -{ P (t ) + }{<, + ^2 y^ 

** 43 y? +♦ 44 y? +^5 ^9^ - 2 u r {(I y^ 


♦ 



^2 ■'■’^ 43 ^6 44 ^8 45 ^10 > “2c y^eJJ; y^ 

d y (B - 9) 

= 'IT? = ^9 (f , y^, , y^) =yio (B - 10) 

d 70 , r 

“ “dT ~ ^10 » yi» ym^ “ Y y^o 

- { P ( T ) + u^} { y ^ + (>^2 73 ■*“ <^53 75 + 4*51^ 7^ 
+ <>^5 79} - 2 u r{ ’I' 51 72 + ’f' 52 ^4 + ’f' 53 ^6 
+ ’I' 54 ^8 +’^55^10^ ^10 - ^*5^9 (B-11) 


The above ten simultaneous first order linear diffe- 
rential Eqs (D-2 to B-11) are solved by Haming’s method 
(Predictor - corrector metaod) 
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The integrals are computed numerically by Guassian- 
Quadrature 15 point formula ♦ The values of these integrals 
were found approximately same by numerical integration as 
well as by emperical formiilae. 
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